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基于 D-Be 紧凑 型 快 中 子 源 硼 中 子 的 慢 化 整形 研究 
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摘要 ” 硼 中 子 俘获 治疗 是 一 种 细胞 尺度 内 强 刘 向 性 与 高 传 能 线 密 度 的 二 元 放射 疗法 ， 具 有 治疗 周期 短 、 对 周转 
素 组 织 损伤 小 等 优势 ， 是 一 种 具有 广阔 前 景 的 癌症 治疗 方法 。 在 BNCT 装置 中 ， 束 流 整形 组 件 的 作用 是 将 
子 源 产生 的 快 中 子 束 慢 化 至 超 热 中 子 能 区 (0.5eV<E<10keV)， 并 保证 中 子 方向 性 。 本 文 使 用 蒙特 卡 罗 模 拟 
序 GEANT4 以 及 FLUKA 模拟 ?Be(d,n)!oB 反应 中 子 源 产生 及 后 续 中 子 慢 化 。 采 用 1.45 MeV、30mA HIIR 
击 9hm 的 薄 皱 靶 ， 由 该 反应 产生 的 中 子 作 为 源 项 ， 对 BSA 进行 了 可 行 性 方案 设计 。 结 果 表 明 : 采用 45cm 
Ý B 下 3 和 5 cm 厚 的 TiF; 组 合 慢 化 层 ，12 cm 厚 的 Pb 反射 层 ，11 cm 厚 的 ALO; 补充 慢 化 层 ， 以 及 0.1mm 
*] Cd 热 中 子 吸收 层 , 保证 出 口 处 的 y 成 分 、 快 中 子 成 分 . Dean 以 及 Que ff IAEA 所 推荐 的 建议 值 。 
研究 得 出 了 低能 气 束 与 薄 皱 丢 的 中 子 谱 及 BSA 具体 设计 方案 ， 为 D-Be 中 子 源 中 子 的 慢 化 整形 提供 了 数据 
参考 ， 为 D-Be 源 进一步 优化 研究 提供 支持 。 
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Abstract [Background]: Boron Neutron Capture Therapy (BNCT) is a binary radiation therapy with strong targeting 
and high energy transfer line density at the cellular scale.It has the advantages of short treatment cycle and minimal 
damage to surrounding healthy tissues,making it a promising cancer treatment method. [Purpose]: This study aims to 
design beam shaping assembly(BSA)to make the neutron beam of D-Be neutron source suitable for BNCT and ensure 
neutron directionality. [Methods]: This article uses Monte Carlo simulation programs GEANTA and FLUKA to simulate 


the generation of?Be(d,n)'9B reaction neutron sources and subsequent neutron moderation.A feasible scheme design for 
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BSA was carried out using a 1.45 MeV,30 mA deuterium beam to bombard a 9 um thin beryllium target,and set a basis 


BSA model with a cylindrical structure as a whole. [Results]: The results show that using a 45 cm thick BiF3 and 5 cm 


thick TiF; combined slowing layera 12 cm thick Pb reflector layer,an 11 cm thick ALO; supplementary slowing 


layerand a 0.1 mm thick Cd thermal neutron absorption layer,the outlet is ensured to y and fast neutron 


composition,O 


epi/ D i, Popi 


/Dest meets the recommended values of the IAEA. [Conclusions]: This study obtained the 


neutron spectra and BSA specific design scheme of low-energy deuterium beams and thin beryllium targets,providing 


data reference for the slowing shaping of neutrons in D-Be neutron sources and supporting subsequent research on D- 


Be sources. 


Key words 


中 子 俘 获 治疗 (Neutron Capture Therapy, NCT) 


于 1936 年 由 G.Locher 首先 提出 , 他 指出 1B 俘获 


Boron neutron capture therapy,D-Be neutron source,Beam Shaping Assembly,Neutron moderation 


一 个 热 中 子 后 会 产生 一 种 高 能 、 


产物 ,这 种 产物 有 较 强 的 生物 效应 ， 因 此 , 中 子 俘 


短 射 程 的 放射 性 


E 作 为 一 种 新 的 癌症 治疗 方法 。 硼 


获 治疗 法 有 可 
中 子 俘 获 治 疗 ( Boron 
Therapy,BNCTO 具有 精准 肿 


Neutron Capture 


Jerg I PE. 183r 5] Hj] 
短 、 对 周围 健康 组 织 损 伤 小 等 优势 , 是 一 种 具有 ) 
阔 前 景 的 癌症 治疗 方法 。 A 
将 携带 10B 的 靶 向 分 子 药物 注射 到 人 体 ，10B 高 特 


明 中 子 俘 获 治 疗 原 理 为 


异性 地 集聚 在 肿瘤 组 织 内 ， 随 后 使 用 经 慢 化 后 的 


低能 量 超 热 中 子 束 对 肿瘤 部 


组 织 内 发 生 OBY Li 反应 , Ù 
子 两 个 射程 仅 约 10 ym〔( 约 一 个 癌 细 胞 大 小 )、 高 


立 进行 外 照射 ， 肿 瘤 


出 wx 粒子 和 TL A 


传 能 线 密度 的 习 


EAT, KIMY 


癌 细 胞 ， 对 肿瘤 周围 健康 组 织 日 


美国 的 W. 


Sweet && ARIF 1 


反应 堆 热 中 


束 的 硼 中 子 俘获 治疗 装置 ， 首 先进 


行 了 硼 中 子 俘获 治疗 临床 试验 ， 
含 硼 药 物 没 有 较 好 地 实现 在 肿瘤 组 织 的 浓 聚 ， 试 
验 中 的 患者 生存 期 均 小 于 一 年 。1958 年 ，Snyder 
等 人 BI 首次 合成 了 硼 茶 丙 氨 酸 (BPA). 日 本 


Mishima 55 A 4l 
将 BNCT 应 用 


于 1987 年 使 


| 胞 水 平定 点 击 杀 
的 损伤 极 小 串 。 
951 年 研究 基于 核 


由 于 当时 使 用 的 


] BPA 作为 硼 化 合 物 


于 恶性 黑色 素 瘤 , 在 临床 上 实现 了 


BNCT 在 中 枢 神 经 系统 以 外 的 
20 世纪 70 FRES, EZAR H 
ES ACIE. SEIL D TRAC 


2x 
2 


药物 等 ， 在 临床 前 下 


肿瘤 治疗 中 的 应 用 


o 
完 层 出 不 


究 中 也 取得 了 许多 积极 结果 。 


2020 年 3 月， 世界 上 第 一 台 加 速 器 BNCT 设备 和 


全 球 首 个 硼 药物 在 日 本 获 批 。 同 年 5 H, 
fg BNCT 疗法 正式 接受 患者 ， 硼 中 子 俘获 疗法 在 


持 建造 一 座 基于 微 堆 的 


世界 上 首次 正式 进入 临床 应 用 。 
国内 对 BNCT 技术 及 其 发 展 极为 重视 。 中 国 
原子 能 科学 研究 院 周 永 戊 院士 于 2009 年 在 北京 主 
医院 中 子 照 射 器 ， 是 世界 
首座 专 为 BNCT 提供 中 子 源 的 核 装 置 。2014 4E 9 


HAE 


月 该 装置 实施 了 我 国 首 例 黑色 素 瘤 患者 的 BNCT 
临床 试 治 , 获得 了 良好 的 治疗 效果 口 。 中 硼 医 疗 联 


合 厦 门 弘 爱 医院 共同 打造 的 加 速 器 BNCT H 
疗 中 心 已 于 2021 年 8 月 完成 安装 并 成 功 打 靶 出 


R, 


动 临 床 进程 。 


并 于 当月 开启 全 性 能 测试 和 动物 实验 , 全 力 推 
2024 年 4 月 睿 斯 科 医疗 与 晋江 市 展 


开 合作 ， 开 展 科 洽 会 展示 了 最 新 研发 出 有 全 自主 
知识 产权 的 基于 加 速 器 的 硼 中 子 俘获 治疗 


( 


accelerator-based Boron 


Neutron Capture 


Therapy, AB-BNCTO 核 系 统 及 精准 放射 治疗 系统 
等 创新 技术 ， 加 速 我 国医 疗 设 备 国产 化 。 


在 危险 性 大 及 安全 隐患 等 因素 ， 限 制 了 反应 堆 
BNCT 的 临床 应 用 ， 反 应 堆 BNCT 中 子 源 逐 渐 退 


由 于 反应 堆 存 在 占 地 面积 大 , 运 维 成 本 高 、 潜 


出 历史 舞台 向 。 随 着 加 速 器 技术 的 进步 , 可 以 通过 
射频 直线 加 速 、 回 旋 加 速 和 高 压 加 速 等 方式 提供 


安全 可 靠 、 使 用 灵活 方便 的 中 子 。 通过 研究 加 速 器 


加 速 粒 子 束 用 击 贰 材料 的 反应 有 : 


"Li(p,n)'Be. 


PC(p,n)'4N, ?T(d,nHe. ?Be(djn)'?^B 等 。 国 内 外 多 


个 组 织 机 构 都 在 推动 AB-BNCT H5 


发 及 临床 应 


用 [81。 


比 ， 


在 稳定 产生 中 子 方面 ， 与 Lip, Be 反应 相 
?Be(d,n)!^B 具有 更 大 的 优势 。 金 属 饶 具有 合适 


的 热 性 能 和 机 械 性 能 ， 争 的 熔点 和 热 导 率 分 别 为 


128 


7°C 和 190W/m-K, 而 金属 锂 的 熔点 和 热 导 率 分 


别 为 180°C 和 84.7W/m-K。 由 7Li(pm7Be 反应 产 
生 的 中 子 作为 源 项 ， 需 要 考虑 到 锂 靶 的 冷却 问题 、 


剩余 核 Be 


(半衰期 为 53 KO. 以 及 锂 台 的 辐 照 损 


fni. K 


] ?Be(d,n) ^B 作为 源 项 ， 皱 具有 高 熔 


点 和 高 热 导 率 ， 有 效 的 水 冷却 可 以 带 走 这 些 热量 ， 


并 且 10B 是 稳定 同位 素 ， 没 有 


料 ， 


余 的 放射 性 。 
本 论文 考虑 争 蔡 代 锂 作为 紧凑 型 中 子 源 靶 材 
并 分 析 低 能 气 束 麦 击 锌 台 后 中 子 的 产 额 、 能 谱 ， 


2 


H 一 种 基于 «Be(d,n)!^B 反应 的 AB-BNCT 装置 


的 BSA 可 行 性 设计 方案 。 
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EIE: 基于 D-Be 紧凑 型 快 中 子 源 硼 中 子 的 慢 化 整形 研究 
1 ”材料 与 方法 FLUKA (FLUktuierende KAskade) 是 一 个 通 
用 的 蒙 卡 粒子 输 运 工具 ， 可 以 模拟 约 60 种 不 同 的 
EESDEBON OR 粒子 ， 包 括 从 1 keV 到 几 千 TeV 的 光子 和 电子 、 
20 世纪 末 ， 有 学 者 提出 了 基于 Beld) B 的 任意 能 量 的 介子 、20 TeV 以 下 的 强 子 以 及 对 应 
中 子 源 ， 基 于 AB-BNCT 产生 治疗 用 中 子 束 外 ,其 的 反 粒 子 、 从 热能 区 往 上 的 中 子 和 重 离子 等 。 应 用 
反应 式 为 : 范围 包括 : 剂量 学 、 探 测 器 设计 、 字 宙 射 线 、 放 射 
9Be+2D 一 1I0B+n+Q (1) 治疗 等 。 本 论文 使 用 FLUKA 对 整个 束 流 整 形体 
IP: Q=4.36MeV。 在 较 低 的 气 束 能 量 下 ， (Beam Shaping Assembly, BSA ) 进行 设计 和 优化 。 
该 反应 也 有 较 高 的 中 子 产 额 ， 并 且 产 生 的 中 子 能 E 
量 较 其 他 中 子 源 更 低 C96 MeVO. 1.45 MeV 的 气 2 ”中 子 源 产生 
AO 9 pm WWE, CERUHHIUAZU 504keV 的 GEANTA ( GEometry ANdTracking) 是 由 欧洲 
能 量 10], 由 于 Q 下 为 4.36 MeV, 剩余 核 10B 的 第 核子 中 心 (CERN) 于 1994 年 开发 ， 是 一 个 以 CH 
6 能 级 ( 表 1) 在 约 1 MeV Wahe FAT 为 底层 语言 , AORERE, 主要 有 以 下 四 种 反应 : 
其 他 核 被 填充 ， 反 应 中 释放 的 大 部 分 能 量 被 消耗 复合 核反应 、 所 核 前 裂 反 应 、 破 裂 反应 以 及 敲 出 反 
在 OB 的 激发 中 。 通 过 采用 注 靶 将 所 有 反应 抑制 在 pum, agg PAKEA 9 um, rh T eicit 
HIC EL UA TF, KÄ AAAA 子 能 被 消除 ， 而 剩 下 3 : E 前 没有 完善 的 评价 核 数据 库 对 该 反应 进行 
的 几 百 keV 就 作为 发 射 中 子 的 动能 ， 导致 发 射 的 研究 ， 需 要 使 用 GEANT4 中 相应 的 物理 模型 计算 
中 子 能 量 非 常 小 00。 因 此 , 采用 薄 靶 增加 了 低能 中 子 产 额 ， 参 考 M.E.Capoulat 等 人 03 的 文献 中 的 
子 产 额 ， 软 化 中 子 谱 ， 降 低 了 后 续 慢 化 的 难度 。 数据 ， 其 中 子 产 额 为 162xl0tnmC-L L45MeV 
" MRR TERARI HT RE, nÉ 1 所 示 ， 其 
jul pisces d im of the 产 额 为 1.82x10" » mC'. 结果 表明 ， 计 算 结果 整 
accessible states of "B 体 上 符合 参考 数据 。 
一 一 -一 一 一 获取 该 中 子 谱 中 子 动量 、 位 置 、 发 射 方向 等 信 
LU Nr A NE 息 ， 作 为 后 续 蒙特 卡 罗 程 序 FLUKA 源 文件 的 输 
Energy level Energy level energy Threshold energy 入 。 
基态 
1 718.380 32 
2 1740.05 E n 
3 2154.27 z 
4 3587.13 pd 
5 4774.0 504.9 220 
6 5110.3 916.4 E T 
7 5163.9 982.0 i 
8 5180 1004.2 s E 
Z 
0 


1.2 蒙特 卡 罗 模 拟 方法 


GEANT4 (GEometry ANdTracking) 是 由 


核子 中 心 CCERN) 于 1994 年 开发 , 是 一 个 以 CH+ 


为 底层 语言 
在 材料 中 的 
样 的 微观 粒 
离子 以 及 


Inf 


E. IMER, FADE 
用 GEANTA 中 相应 的 物理 


， 以 蒙特 


输 运 过 程 的 工具 ， 


FE 罗 方法 为 基础 的 模拟 粒子 


AA 


t B6 


f ORTA P3 


可 衰变 核 素 的 
MP AS AST 


AEN 


欧洲 


够 涉及 各 种 各 


、 质 子 、X 射线 、 各 种 
辐射 粒子 )， 主 要 在 核 物 
页 域 中 使 用 。 本 论文 使 
模型 模拟 D-Be 反应 。 


H 


图 1 1 


全 
in 


.0 
0 lxl056 2x10 3x10 4x10 5x10  6x105 

Neutron energy/eV 

.45 MeV ATO kil UBI. 9Be(d,n)!?B 反应 的 中 


子 谱 


Fig.l Neutron spectrum of simulated ?Be(d,n)'9B reaction 


with 1.45 MeV deuterium beam bombard on a thin target 
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BNCT 装置 设计 
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束 流 整形 体 BSA 


工作 的 核心 内 容 之 一 是 对 


E 技 Ñ 
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BSA 进行 设计 和 优化 ， 其 作用 是 将 快 中 子 慢 化 至 
可 用 于 治疗 的 固定 能 量 范 围 的 超 热 中 子 。 在 BSA 
优化 中 ， 每 个 组 件 都 应 该 根据 上 游 组 件 的 中 子 输 
出 进行 优化 。 经 慢 化 后 BSA 出 口 处 中 子 束 流 的 相 
关 参 数 应 满足 国际 原子 能 机 构 C Intemational 
Atomic Energy Agency, IAEA ) 所 给 出 的 建议 值 04， 
如 表 2 所 示 。 


表 2 IAEA-TECDOC-1223 推荐 的 治疗 用 中 子 束 指标 
Table 2 Targeted values of therapeutic neutron beam 
recommended by IAEA-TECDOC-1223 


BNCT 端口 参数 建议 值 

BNCT port parameters Recommended values 
Doi (nccm?:s) -10? 

D,. /® pi (Gy:cm?) <2x10-13 

D, /Dy (Gy:cm?) <2x10-3 

Dopi/ Dest >20 

OD, >20 

J/® ~0.7 


快 中 子 能 区  E/eV 

Fast energy group 

超 热 中 子 能 区 E/eV 
Epithermal energy group 
热 中 子 能 区 —E/eV 


Thermal energy group 


E>10000 


0.5«E«10000 


E«0.5 


其 中 ， 快 中 子 成 分 Ds /Pepi M y BG) D, /De 
分 别 为 快 中 子 和 y 的 吸收 剂量 率 与 超 热 中 子 通 量 
的 比值 ; Deygu 为 超 热 中 子 与 快 中 子 通 量 的 比值 ; 
Ooi/Ons 为 超 热 中 子 与 快 中 子 通 量 的 比值 ; 
@Desy Dtotal 为 快 中 子 与 总 中 子 通 量 的 比值 ， 流 量 通 
量 比值 J@ 表征 了 中 子 束 的 前 向 1 


3.1 BSA 结构 


设 定 一 个 基准 BSA 模型 ,整体 为 圆柱 体 结构 ， 
如 图 2 所 示 ， 从 左 至 右 依次 为 : 反射 层 、 气 束 通 
道 、 靶 体 、 冷 却 体 、 慢 化 层 、 准 直 体 、 热 中 子 吸收 
层 。 气 束 从 左 侧 注 入 , 靶 材 料 外 层 沉积 在 相对 较 薄 
的 层 〈 由 钒 、 锟 等 支撑 材料 ) 上 ， 以 避免 氧 脆 。 本 
论文 中 不 考虑 这 些 背 衬 材料 ， 因 为 它们 对 模拟 的 
结果 几乎 没有 影响 。 皱 靶 后 端 为 铜 腔 ， 厚 度 为 1.5 
cm， 铜 腔 中 舍 水 ， 有 效 的 水 冷却 能 够 安全 地 带 走 
束 流 沉积 的 热量 05]。 准 直 孔 的 出 口 半径 为 7.5 cm. 
根据 D-Be 反应 中 子 谱 具 体 信 息 及 基准 模型 , 使 用 
FLUKA 对 整个 BSA 进行 可 行 性 方案 设计 。 


HE 


NA! 

Beam channel 
Bi Target 
H Copper chamber 
E water 
Bl Thermal neutron 
absorption layer 


Reflector layer 


图 2 基准 BSA 模型 


Fig.2 Basis BSA model 


3.0 $KE 


慢 化 层 的 主要 作用 是 将 快 中 子 尽 可 能 多 地 慢 
化 成 超 热 中 子 。 慢 化 层 材 料 应 在 超 热 中 子 区 具有 
低 散 射 截面 , 在 快 中 子 区 具有 高 散射 截面 , 并 对 超 
热 中 子 的 吸收 截面 应 尽 可 能 低 09。 本 论文 采用 多 
层 慢 化 组 合 结构 。 由 于 含 氟 材 料 具 备 良 好 的 慢 化 
快 中 子 性 能 ， 选 择 PbpF2 等 含 重 元 素 的 材料 作为 第 
一 层 慢 化 材料 ，MgF> 等 含 较 轻 元 素 的 材料 作为 第 
二 层 慢 化 材料 [7]。 

参考 众多 学 者 016-19 对 慢 化 组 合 层 的 选择 及 模 
拟 , 将 慢 化 层 组 合 结构 总 长 设置 为 50 cm。 第 一 层 
慢 化 层 厚 度 变 化 区 间 为 15-45 em; 步 长 设 为 Scm， 
半径 为 25 cm。 通过 对 不 同 厚度 慢 化 材料 组 合 进行 
模拟 , 得 到 候选 材料 与 厚度 的 最 佳 组 合 。 计 算 结果 
如 图 3、 图 4、 图 5 所 示 ， 分 别 给 出 了 不 同 组 合 下 
的 中 子 慢 化 效果 。 


着 


一 一 BiF;+AIF。 一 * 一 PbF,+AIF， 一 * 一 PbFI+AIF， 
一 * 一 BiF;+MgF， 一 = PbF,-MgF, —— PbF,+MgF, 
一 一 BiF;+TiF。 一 * 一 PbF,+TiF。 —e— PbF,+TiF, 


Thickness of the first slowing layer/cm 


图 3 不 同 慢 化 层 组 合 下 040, B8 EE EAE (S 
Fig.3 Relationship between the variation of 0/0, with 


thickness under different combinations of slowing layers 
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$, /10" mcm? s! 


FA BiF;FAIF, —9— PbF,+AlF, —4— PbF,+AlF, 
3,0 |—e— BiF;*MgF, —— PbF,+MgF, —A— PbF,*MgF, 
C——^— BiF,*TiF, —sm— PbF,+TiF, —e— PbF,+TiF, 


15 20 25 30 35 40 45 
Thickness of the first slowing layer/cm 
图 4 不 同 慢 化 层 组 合 下 @si 随 厚度 变化 关系 
Fig.4 Relationship between the variation of ® 


epi with 


thickness under different combinations of slowing layers 


—— BiF,-AIF, —5— PbF,-AIF, —4— PbF,+AlF; 
—e— BiF,+MgF, —5— PbF,4MgF, —— PbF4MgF, 
— — BiF,4TiF, —s5— PbF,+TiF, —e— PbF,4TiF, 


15 20 25 30 35 40 45 
Thickness of the first slowing layer/cm 
图 5 不 同 慢 化 层 组 合 下 @aswVdiom 随 厚度 变化 关系 
Fig.5 Relationship between the variation of Des/ Diota with 


thickness under different combinations of slowing layers 


由 图 3 和 图 5 可知, BiFs 和 MgFo 的 组 合 层 的 
dey@ast 较 其 他 组 合 层 高 , Drosy/ Dota 较 其 他 组 合 层 
低 。 第 一 层 慢 化 层 厚 度 为 45 cm 时 ，BiF3 和 TiF3 
HEE @ey@ast 比值 达 最 大 , Dasy Dwotal 达 最 低 , 并 
由 图 4 nA ei 随 BiFs 厚度 的 增加 而 增 大 。 选 
择 25 cm 厚 的 BiFs 和 25 cm 厚 的 MgF2 与 45 cm 
E E BiF3 和 5 cm 厚 的 TiFs 作为 候选 慢 化 层 组 合 

案 。 中 子 经 过 这 两 种 慢 化 层 组 合 后 的 能 谱 如 图 6 
pod 


—— BiF,«TiF, 
7 — BiFj4MgF, 


6 nem? s! 
Pota 10° nr ems 
w A Un a 


N 


0 - das 
10? 10! 10! 10? 10 10” 
Neutron energy/eV 
图 6 中 子 经 过 两 种 慢 化 层 组 合 后 的 能 谱 
Fig.6 Energy spectrum of neutrons after passing through 


two types of slowing layers 


由 图 6 可 知 ， 采 用 厚度 为 4 em 的 BiF3 和 5 
em 的 TiFs 慢 化 组 合 层 , 超 热 中 子 范围 内 的 中 子 更 
加 集中 ， 热 中 子 通 量 更 低 。 厚 度 为 25 cm 的 BiFs 
和 25 cm 的 MgFo 组 合 层 结构 ， 其 热 中 子 通 量 较 
高 。BSA 设计 目标 其 一 就 是 得 到 很 低 的 出 口 处 热 
中 子 通 量 ， 因 此 ， 采 用 厚度 为 45cm 的 BiF3 和 5 
cm 的 TiFs 作为 慢 化 组 合 层 。 


3.3 反射 层 


反射 层 的 主要 作用 减少 中 子 泄漏 ， 最 大 化 
BSA 出 口 处 的 中 子 通 量 。 反射 体 为 圆 简 状 ， £8 
在 中 子 源 和 慢 化 体外 侧 ， 采 用 在 超 热 中 子 能 区 具 
有 高 弹性 散射 截面 和 低 吸 收 截 面 的 铅 、 锐 两 种 反 
射 层 材料 , 针对 不 同 厚 度 进 行 模 拟 。 计 算 结 果 如 图 
7 和 图 8 所 示 。 


Lead 1 
Bi q 


Reflector thickness/cn 
图 7 铝 、 锐 反射 层 下 Guess 随 厚 度 变化 关系 
Fig.7 Relationship between the variation of Doi/Dest with 


thickness under lead and bismuth 
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o, ,/107 ncm? s? 


0 5 10 15 20 25 30 


Reflector thickness/cm 
图 8 £i. BRIEF @。 随 厚度 变化 关系 
Fig.8 Relationship between the variation of D, with 


thickness under lead and bismuth 


由 图 7 和 图 8 可 知 , 增加 一 定 的 反射 层 厚 度 ， 
能 够 使 得 Depi/Drast LEIE, Depi 明显 增 大 ， 且 铅 的 效 
果 要 优 于 馈 。 本 论文 研究 适合 安装 在 医院 的 紧凑 
型 快 中 子 源 装 置 , 适当 的 选择 反射 层 厚度 , 可 以 缩 
小 整体 装置 的 尺寸 。 因 此 ， 本 论文 选择 厚度 为 12 
em 的 铅 作 为 反射 层 。 


3.4 补充 慢 化 层 及 准 直 器 


由 于 经 过 慢 化 层 组 合 结构 的 中 子 并 未 完全 慢 
化 至 所 需 超 热 中 子 能 量 范围 ， 需 要 通过 补充 慢 化 
层 进一步 慢 化 。 采 用 含 较 轻 元 素 的 AbOs. AIFs. 
MgO 这 3 种 材料 针对 不 同 的 厚度 进行 模拟 计算 。 
补充 慢 化 层 为 圆柱 结构 , 位 于 慢 化 层 组 合 下 端 , 半 
径 与 慢 化 层 组 合 保持 一 致 ， 其 厚度 变化 区 间 为 5- 
23 cm， 步 长 设 为 3cm， 计 算 结 果 如 图 9、 图 10 所 
不 。 


Thickness/cn 
图 9 不 同 慢 化 层 下 Quo Bt EE EAE IOS RR 
Fig.9 Relationship between the variation of Doi/Dest with 


thickness under different slowing layers 
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9 —4— MgO | ] 
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Thickness/cm 
图 10 不 同 慢 化 层 下 Dw; 随 厚度 变化 关系 
Fig.10 Relationship between the variation of D, with 


thickness under different slowing layers 


由 图 9 可 看 出 ， 随 着 三 种 补充 慢 化 层 材料 厚 
度 的 增加 ，@epy@ast 比值 也 随 之 增 大 ， 厚 度 为 14 
cm 左右 时 ,三 种 材料 的 @enyg@as 比值 大 于 20。 由 
图 10 可 以 看 出 ， 随 着 三 种 补充 慢 化 层 材 料 厚 度 的 
增加 ， 超 热 中 子 通 量 Doi HEZ R. BIA MgO M 
料 能 够 在 更 低 的 厚度 下 达到 @eoygGast 比值 要 求 ， 
但 是 其 @ei 值 始终 比 AIF 的 Boni 值 低 。 根 据 BSA 
设计 要 求 , 在 满足 @ey@ast 要 求 的 同时 , 选择 超 热 
中 子 通 量 更 大 的 补充 慢 化 层 材料 AF, REEN 
17 em. 
因 束 流 整 形体 BSA 出 口 处 的 中 子 较为 发 散 ， 
准 直 体 一 般 采 用 锥 形 孔 道 结构 以 引导 中 子 束 流 至 
患者 肿瘤 部 位 进行 照射 治疗 。 将 厚度 为 17 cm 的 
AIF; 代入 准 直 器 锥 形 结构 中 ， 其 内 径 为 25 em, 
外 径 设 置 为 7.5 cm。 计 算 结果 : Ouen 比值 为 
27.18. EH 10 可 知 ,， 可 以 设置 更 小 厚度 的 补充 慢 
化 层 材料 以 增 大 超 热 中 子 通 量 。 因此 , 选择 补充 慢 
化 层 材料 AlFs， 其 厚度 为 11 cm。 计 算 结 果 : 
Dop/ Drast 比值 为 21.79， 满 足 要 求 。 


3.5 热 中 子 吸 收 层 


经 过 上 述 慢 化 层 组 合 及 补充 慢 化 层 的 方案 下 ， 
不 可 避免 的 会 产生 热 中 子 ，Dop/@u 比值 为 9.19， 
不 满足 IAEA 最 小 限 值 20 的 要 求 。 热 中 子 会 对 人 
体 皮 肤 和 表层 组 织造 成 伤害 ， 因 此 在 超 热 中 子 束 
照射 之 前 ， 需 添加 一 层 热 中 子 吸收 层 以 降低 热 中 
子 比例 。 该 材料 应 对 能 量 小 于 0.5eV WARTE 
有 很 大 的 吸收 截面 ， 在 中 能 区 和 高 能 区 具有 很 低 
的 散射 截面 和 吸收 茂 面 。 

常见 的 热 中 子 吸收 材料 为 Li 和 Cd, SLi 对 
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中 子 的 吸收 截面 为 940b， 而 且 Li(n,W) 反 应 吸收 
中 子 产 生 的 Y 污染 很 少 ,不 会 造成 后 续 处 理 难 度 。 
Scd 对 中 子 的 吸收 截面 为 2100b， 中 能 区 平均 热 
中 子 吸 收 截 面 为 8b， 高 能 区 平均 热 中 子 吸收 截面 
为 0.05 bp， 其 自然 丰 度 为 12.2%。 

采用 LiF、 天 然 锅 这 两 种 材料 针对 不 同 的 厚度 
进行 模拟 计算 。 热 中 子 吸收 层 为 圆 饼 状 薄 层 结 构 ， 
紧 靠 补充 慢 化 层 ， 热 中 子 吸收 层 厚 度 变化 区 间 为 
0.01-0.15 cm, 步 长 设 为 0.01 em, 计算 结果 如 图 11 
所 示 。 


40 


0.06 


0.08 — 0.10 


Thickness/cm 
图 11 不 同 热 中 子 吸 收 层 下 Dwi/DPw 随 厚度 的 变化 关系 


Fig.11 Relationship between the variation of $,,/0, with 


thickness under different thermal neutron absorption layers 


由 图 11 可 以 看 出 , LIF 材料 的 duda 比值 随 
热 中 子 吸收 层 厚度 的 增加 而 缓慢 增加 ， 变 化 幅度 
较 小 。LiF 厚度 为 0.07cm Hf, Depi/On 为 20.4。Cd 
材料 的 edu 比值 随 热 中 子 吸收 层 厚度 的 增加 而 
剧烈 增加 , 初始 厚度 为 0.01 cm 时 Depi/ On X 39.03, 
厚度 为 0.14 cm 时 dea 为 610.9， 升 高 了 约 
1482%, 吸收 热 中 子 效果 十 分 显著 。 因 此 , 在 本 论 
文中 , 使 用 0.01 cm 的 锣 热 中 子 吸收 层 就 可 达到 要 
R, Depi/Dn 比值 为 39.03。 此 时 , BSA 出 口 处 中 子 
RAHI y RIN D, /Oepi JI 1.24x1075 Gy'cm2， 快 中 
子 成 分 DD,, /Doi 为 3.99x10-14 Gy-cm?， 小 于 IAEA 
最 大 限 值 2x10-3 Gy:cm?， 满 足 要 求 。 


4 ”结论 
AJRAYE OC HIORGSGH SEU, H ?Be(d,n)'?B 


反应 产生 的 中 子 作为 源 项 , 产 额 达 1.82x10!! ncmC- 
1 经 对 束 流 整 形体 进行 设计 和 优化 可 得 BSA 优化 
设计 方案 : 采用 45 cm 厚 的 BiF3 和 5 cm EH TiFs 
组 合 慢 化 层 ，12 cm EH Pb 反射 层 ，11 cm 厚 的 
ALbOs 补充 慢 化 层 , 以 及 0.1 mm 厚 的 Cd 热 中 子 吸 
收 层 ,空气 端 出 口 参数 Oei 7J 1.05x 108 ncm? s1, 
Dop/Dt 为 39.03, Pepi/Drast 为 21.82，D, /Teri 为 
1.24x10-15 Gy:em?, D, /Dopi 为 3.99x10-14Gy:cm2， 
J/@ 约 为 0.6。 

经 该 设计 方案 下 BSA 整形 慢 化 后 的 超 热 中 子 
束 流 能 谱 如 图 12 所 示 ,， 超 热 中 子 束 流 的 参数 值 和 
IAEA 建议 值 以 及 世界 上 其 它 超 热 中 子 BSA 设计 
的 对 比如 表 3 所 示 。 在 保证 出 口 处 的 y 成 分 、 快 
FP BEAT. Depi/On 以 及 Pepi/Orast 符合 IAEA 所 推 
荐 的 控制 值 的 前 提 下 ， 当 电流 大 于 0.28A 可 满足 
dei>1x10?n'cm-2.s-1。 本 研究 为 D-Be 中 子 源 中 子 
的 慢 化 整形 提供 了 数据 参考 , 为 D-Be 源 的 进一步 
优化 研究 提供 支持 。 

1.5 
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Neutron energy/eV 
图 12 经 慢 化 后 的 超 热 中 子 束 流 能 谱 
Fig.12  Epithermal neutron beam energy spectrum after 


moderation 


表 3 最 终 BSA 出 口 处 中 子 束 流 相 关 参 数 计算 值 和 其 它 设计 方案 对 比 


Table 3 Comparison between the calculated values of neutron beam parameters at the final BSA outlet and other design 


schemes 

BNCT 端口 参 关 D(X109 D pa /0,,, D, /$,.. 

T] 数 epi Dan/ Dp De/ Dn fast / * epi y ! * epi Jio 
BNCT port parameters n'cm? s!) (x10:8 Gy'cm2) (X10? Gy:cm?) 
建议 值 

~l >20 >20 «2 «2 ~0.7 

Recommended value 
本 论文 This work 0.11 21.82 39.03 0.0124 0.399 0.6 
Cerullo 等 09] 0.66 23.2 133 3.19 1.1 0.58 
Montagnini 等 [20] 1.23 35.0 103.1 1.7 1.7 0.61 
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表 3 最 终 BSA 出 口 处 中 子 束 流 相 关 参 数 计算 值 和 其 它 设计 方案 对 比 CHE) 


Table 3 Comparison between the calculated values of neutron beam parameters at the final BSA outlet and other design 
schemes(Continued table) 


BNCT 端口 参数 Deni(X109 人 
BNCT port parameters ”ncm2.S-D epi fost 
Fantidis 等 20 1.096 - 

ZHU Yinan 等 (7 1.91 21.1 


epi 


D, //6,. D, Dai 
JO fast epl YN: epl JỌ 
*^ (x10-3 Gy:cm?) (X103 Gy:cm2) j 
178.6 1.4 1.4 z 
114 0.99 1.79 0.96 


作者 贡献 声明 ” 王 佳 负责 模型 构建 与 仿真 建 模 、 
文章 起 草 与 修订 ; 刘 义 保 负责 研究 设计 、 理 论 指 导 、 


审阅 文章 内 容 与 修订 。 魏 强 林 负责 研究 设计 、 理论 


已 


Hus Kd HE. FIE EEREN RN. 


参考 文献 


LUDEWIGT B A, BLEUEL D L. Clinical requirements 
and accelerator concepts for BNCT[R]. Accelerator and 
Fusion Research Division,LBNL-40262Rev UC- 
406,2,US:Lawrence Berkeley National Laboratory. 
1997. 

FARR L E, SWEET W H, LOCKSLEY H B, et al. 
Neutron capture therapy of gliomas using boron[]J]. 
Transactions ofthe American Neurological Association, 
1954, 13(79th Meeting): 110-113. PMID: 13238328. 
SNYDER H R, REEDY A J, LENNARZ Wm J. 
Synthesis of Aromatic Boronic Acids. Aldehydo 
Boronic Acids and a Boronic Acid Analog of Tyrosine 
1[J/OL]. Journal of the American Chemical Society, 
1958, 80(4): 835-838. DOI:10.1021/ja01537a021. 
MISHIMA Y, ICHIHASHI M, HATTA S, et al. First 
human clinical trial of melanoma neutron capture. 
Diagnosis and therapy[J]. Strahlentherapie Und 
Onkologie: Organ Der Deutschen 
Rontgengesellschaft ... [et AI], 1989, 165(2-3): 251-254. 
PMID: 2494743. 

YONG Z, SONG Z, ZHOU Y, et al. Boron neutron 
capture therapy for malignant melanoma: first clinical 
case report in China[J/OL]. Chinese Journal of Cancer 
Research, 2016, 28(6): 634-640. 
DOI:10.21147/j.15sn.1000-9604.2016.06.10. 

MENG P, HE G Z, ZHANG Q W, etal. Study ofneutron 
production and moderation for sulfur neutron capture 
therapy [J/OL]. Nuclear Science and Techniques, 2019, 
30. DOL:10.1007/s41365-018-0529-3. 

KIYANAGI Y, SAKURAI Y, KUMADA H, etal. Status 
of accelerator-based BNCT projects 
worldwide[C/OL /AIP Conference Proceedings. AIP 


Conf. Proc. 2160, 050012 (2019). 

8 ALTIERI S, PROTTI N. A brief review on reactor- 
based neutron sources for boron neutron capture 
therapy[J/OL]. Therapeutic Radiology and Oncology, 
2018, 2: 47-47. DOL10.21037/tro.2018.10.08. 

9 COLONNA N, BEAULIEU L, PHAIR L, et al. 
Measurements of low-energy (d,n) reactions for BNCT. 
Boron Neutron Capture Therapy[J/OL]. Medical 
Physics, 1999, 26(5): 793-798. DOI:10.1118/1.598599. 

10 | CAPOULAT M E, HERRERA M S, MINSKY D M, et 
al — 9Be(d,n)lOB-based neutron sources for 
BNCT[J/OL]. Applied Radiation and Isotopes, 2014, 88: 
190-194. DOL:10.1016/j.apradiso.2013.11.037. 

11 CAPOULAT M E, SAUZET N, VALDA A A, et al. 
Neutron spectrometry of the 9Be(d (1.45 MeV), n)10B 
reaction for accelerator-based BNCT[J/OL]. Nuclear 
Instruments and Methods in Physics Research Section 
B: Beam Interactions with Materials and Atoms, 2019, 
445: 57-62. DOL:10.1016/j.nimb.2019.03.005. 

12 ” 张 雅 玲 . 9Be(d,n) fll 9Be(pyn) 反 应 加 速 器 中 子 源 中 子 
学 特性 研究 [D/OL]， 兰 州 : 兰州 大 学 , 2018. 


Zhang Yaling. Research on neutronics characteristics of 


9Be(dn) and 9Be(p,n) reaction accelerator neutron 
sources[D/OL]. Lanzhou: Lanzhou University, 2018. 

13 | CAPOULAT M E, MINSKY D M, KREINER AJ. 
Computational assessment of deep-seated tumor 
treatment capability of the 9Be(d,n)1 OB reaction for 
accelerator-based Boron Neutron Capture Therapy 
(AB-BNCT)[J/OL]. Physica Medica, 2014, 30(2): 133- 
146. DOL:10.1016/j.ejmp.2013.07.001. 

14 AGENCY I A E. Current Status of Neutron Capture 
Therapy: IAEA-TECDOC-1223[R]//Current Status of 

Neutron Capture Therapy. International Atomic Energy 

Agency, 2001: 29-31. 

15 NAKAYAMA S. Simulated performance evaluation of 

d-Be compact fast neutron source[J/OL]. Journal of 

Nuclear Science and Technology, 2023: 1-7. 

DOI:10.1080/00223131.2023.2225498. 

16 ” 王 勇 泉 . 加 速 器 硼 中 子 俘获 治疗 装置 束 流 整 形体 的 


XXXXXX-8 


王 佳 : 基于 D-Be 紧凑 型 


INS 


FP 子 源 硼 中 子 的 慢 化 整形 研究 


17 


18 


19 


20 


21 


设计 及 其 临床 剂量 学 参数 研究 [D/OL]. 兰州 : 兰州 
大 学 , 2023. 
Wang Yongquan. Design of Beam Shaping Assembly 


for accelerator boron neutron capture therapy device 
and study of clinical dosimetric parameters [D/OL]. 
Lanzhou: Lanzhou University, 2023. 

朱 益 楠 ， 林 作 康 ， 户 林 远 ,等 . 基于 和 气 气 中 子 源 硼 中 
子 爷 获 治疗 的 中 子 慢 化 整形 研究 [J]， 核 技术 , 2022, 
45(1): 29-36. DOIL10.11889/.0253- 
3219.2022.hjs.45.010202. 

Zhu Y-N, Lin Z-K, Lu L-Y, et al. Design of beam 


shaping assembly for boron neutron capture therapy 
basedon D-T neutron source[J]. Journal of Nuclear 
Science and Technology, 2022, 45(1) 29-36. 
DOIL:10.11889/j.0253-3219.2022.hjs.45.010202. 
HAIR, 胡 志 良 ， 童 剑 飞 ， 等 . 基于 3.5 MeV 射频 
极 质子 加 速 器 硼 中 子 俘 获 治疗 装置 的 束 流 整形 体 设 
计 [J/OL]. 物理 学 报 , 2019, 67(14): 142801-142801. 
DOI:10.7498/aps.67.20180380. 

Tian Y-S, Hu Z-L, Tong J-F, et al. Design of beam 


zd 


shaping assembly based on 3.5 MeV radio-frequency 
quadrupole proton accelerator for boron neutron 
capture therapy[J/OL]. Acta Physica Sinica, 2019, 
67(14): 142801-142801. 
DOI:10.7498/aps.67.20180380. 

CERULLO N, ESPOSITO J, LEUNG K N, et al. An 
irradiation facility for Boron Neutron Capture Therapy 
application based on a radio frequency driven D-T 
neutron source and a new beam shaping 
assembly[J/OL]. Review of Scientific Instruments, 
2002, 73(10): 3614-3618. DOI:10.1063/1.1505128. 

MONTAGNINI B, CERULLO N, ESPOSITO J, etal. 
Spectrum Shaping of accelerator based neutron beams 
for BNCT[J]. Nuclear Instruments and Methods in 
Physics Research A, 2002, 476: 90. 
DOI:10.1016/S0168-9002(01)01398-5. 

FANTIDIS J G. Beam shaping assembly study for 
BNCT facility based on a 2.5 MeV proton accelerator 
on Li target[J/OL]. Journal of Theoretical and Applied 
Physics, 2018, 12. DOI:10.1007/s40094-018-0312-1. 


XXXXXX-9 


